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^J ' Abstract 

The Casimir energy of quantum fluctuations about the classical kink configuration is com- 
puted numerically in the weak coupling approximation for a recently proposed lattice sine- 
Gordon model. This energy depends periodically on the kink position and is found to be 
approximately sinusoidal. 
PACS classification numbers: 03.65.Sq, ll.10.Lm, 63.10.+a. 
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^ ! 1 Introduction 

[ — ' Classically, a topologically stable field configuration may be regarded as lying in a potential well in 

OS ■ an infinite dimensional configuration space. Two solutions in different sectors are separated by an 

infinite potential barrier by virtue of the topology. Quantum mechanically, a particle cannot sit at 

the bottom of a well: it always possesses a zero point energy dependent on the shape of the well 

Q_i' bottom. The analogous situation applies to fields also. 

We wish to find the zero point energy associated with a kink configuration in a certain lattice 
sine-Gordon model. The model was chosen because it has an exact lattice version of the topological 
lower bound on kink energy, and an explicit static kink which saturates this bound. The calculation 
is performed in the weak coupling approximation, by performing a Taylor expansion of the potential 
about the kink configuration, truncating the expansion at quadratic order. We then find the normal 
modes of the truncated system, reducing the problem to that of an infinite set of decoupled harmonic 
oscillators. On quantization, each contributes a zero point energy. Summing over all oscillators gives 
an infinite total. 

If we similarly approximate the vacuum well-bottom we can derive the zero point energy associ- 
ated with the trivial vacuum. This is also infinite. The "physical" quantity required is the energy 
difference between the kink and the vacuum, since we can always define the vacuum energy to be 
zero. This quantity, analogous to the Casimir energy of Quantum Electrodynamics |l[, turns out to 
be finite and less than the classical energy of the kink. Since the model possesses only discrete trans- 
lation symmetry, the normal mode frequencies, and hence the Casimir energy depend periodically 
on the kink position. 

Calculation of the normal mode frequencies amounts to finding the eigenvalues of an infinite- 
order, tridiagonal, symmetric matrix. In practice this is not possible and the system must be 
truncated symmetrically about the kink centre, ignoring the large \x\ degrees of freedom. It is 
possible to show that the resulting quantum energy correction must be negative, for any size of 
truncated system. The correction may be calculated numerically. 

'Present address: Department of Mathematics, University of Texas at Austin, Austin TX 78712, USA 



2 The Lattice Sine-Gordon Model 

Consider the lattice sine-Gordon model (LSGM) denned by the Lagrangian 
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where /i is the lattice spacing, to is a mass parameter, A is a dimensionless coupling constant and Z 
denotes the set of integers. Taking the h — > limit one recovers the standard sine-Gordon model: 
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where x = hj . The standard sine-Gordon Lagrangian is more usually written in terms of <f> = 1i\) 
and differs by a factor of 8, but the above form is more convenient. 
The LSGM has an infinite set of discrete potential minima, 
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where n S Z, and a static kink interpolating between neighbouring minima, which may be derived 
using a Bogomol'nyi argument. The energy of a static configuration is just the potential V, and 
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when kink boundary conditions are imposed. This (the Bogomol'nyi) bound is saturated if and only 
if 

sin o (^i+i -i>i) = hm sin -(ip j+ i +ipj). (5) 

The first order difference equation, rtq), is called the Bogomol'nyi equation and, remarkably, has 
an explicit kink solution: 
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A 



2 + hm 
2 — hm 



j-b/h 



(6) 



The dimensionless parameter /im £ (0, 2) for sensible solutions. The arbitrary constant b may take 
any real value — the kink energy is not position dependent. The right hand side of (0) may be 
interpreted as the classical kink mass, M . Small velocity kink dynamics can then be approximated 
by geodesic motion of a point particle, mass M, on the manifold generated by b translation, with a 
natural induced metric M. 



3 The Weak Coupling Approximation 

We follow the method outlined in pjadapted for the infinite lattice. We work in natural units, 
effectively having absorbed a factor \jh into A so that, by the standard argument, the weak coupling 
criterion A <C 1 yields a semi-classical approximation. 



The momentum conjugate to ipj is 
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Thus, the LSGM Hamiltonian is 
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Let tp be a static configuration giving a local (in configuration space) minimum of the potential, 
V(tp). We treat motion about this stable configuration in the small A approximation by Taylor 
expansion of V: 

vty) = v$) + ^J2 w ^j - ^•XV'fc -&) + ■■■ (io) 
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Note that W is a real, symmetric, tridiagonal matrix. Note also that the next term in the expansion 

is 
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and that the three derivatives with respect to ip introduce a factor of A 3 , leaving an overall factor of 
A after cancelling the 1/A 2 of V(ip). Thus, higher corrections are at least of order A, which is why 
one can truncate the series in the small A limit. 

Owing to the symmetry of the W matrix, there exists an orthogonal transformation R such that 
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where U is a diagonal matrix. We may reduce the system to a sequence of decoupled harmonic 
oscillators by transforming to the rotated coordinates (normal coordinates) 



£j =yRjk(ipk-ipk), 
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which have conjugate momenta 
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Then, 
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where Q 2 are the eigenvalues of the W matrix, none of which can be negative since ip locally 
minimizes the potential. 

We now quantize in standard canonical fashion, taking ip to be first the vacuum, then the kink 
located at x — b. The vacuum ground state energy is 
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The kink, by virtue of the zero mode b of (0) lies not in a potential well, but in a level-bottomed 
valley meandering through configuration space. One of the normal modes is locally tangential to the 
valley bottom and consequently has vanishing corresponding eigenvalue (zero frequency). We shall 
treat this translation mode, b, classically because in the weak coupling approximation it is much 
heavier than the orthogonal modes (mass to/ A 2 compared with m). 



While the translation orbit of the static kink (0) is an equipotential curve, neighbouring orbits 
are not: the potential varies periodically along them. So the eigenvalues of the W matrix (Qf) 2 will 
be b dependent with period h. The ground state energy of a kink at b is 

E ^ = % + ^kX n ' m - (18) 

3 

where the sum may be taken over all eigenvalues since the one we wish to omit is zero anyway. 

There is no reason to expect either E? ~, (b) or -E? i to converge to a finite sum but one would 
expect finite Casimir energy 

£ = Ef 0} (b)-El Q} -^ (19) 

since the lattice spacing h has effectively cut off the ultra violet divergence problems of the con- 
tinuum model. We must still take care when manipulating the divergent sums of (|19|) that they 
are suitably regulated before being combined. The method of regulation is determined by practical 
considerations: we calculate E? -,(b) and -E? } on a finite lattice, compute £ and then allow the 
lattice size to grow large. In practice we must truncate the lattice to finite size anyway in order to 
solve the kink matrix eigenvalue problem. 

4 The Eigenvalue Problem 

We now address the problem of finding the eigenvalues of the vacuum and kink W matrices, W° 
and W K . Let us first consider the vacuum matrix, obtained from dill) by substituting ipk = 0: 
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This is a very simple matrix. All the diagonal entries are equal, as are all the upper and lower 
diagonal entries. Consequently, the spectrum of the truncated matrix of order N is known exactly 
(see, for example, Q): 

^.o\2 4: + h 2 m 2 4-h 2 m 2 ( jir \ 

<°S> = -^ ^r cos \—i) • (21) 

j = 1, 2, . . . , N. In the N — > oo limit this is strikingly similar to the dispersion relation for phonons 
on the lattice H because the normal modes of oscillation about the vacuum are phonons. 

The question of making sense of the N — > oo limit of J2i ^°j ^ s irrelevant because there is no hope 
of finding the exact spectrum of the kink W matrix, W K . Substituting (g) into (^lj) we obtain an 
explicit expression for W K (b). Since the kink is highly localized (provided hm is not too small), at 
large \j\ the kink configuration rapidly approaches neighbouring vacuum minima so that WK,(b) tend 
to W° k away from the matrix centre. The suggestion, then, is that we truncate the lattice, pinning 
all the large \j\ degrees of freedom, for which the kink and vacuum configurations are essentially 
identical, to their classical values. The resulting finite problem may be solved numerically. 

One can prove that, whatever the order, N (an odd integer), of the symmetrically truncated 
system (that is, truncated symmetrically about the lattice site j = 0), the Casimir energy must be 
negative. The proof rests on the observation that, truncated to order N = 2n + 1, 

W K = W° + B + C, (22) 



where B is the symmetric, tridiagonal matrix, 
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> 0, 
Pj = bj + bj+i, 

and C is the diagonal matrix with elements, 
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Now, C manifestly has negative semi-definite eigenvalues, and B is a special case of a class of matrices 
whose spectra are known to be negative semi-definite ||. Regarding B and C as perturbations to 
the matrix W° , we apply a corollary of the minimax theorem H : 

Theorem: If G,T,T are symmetric N x N matrices with eigenvalues 9 r ,v s ,^t (all sets 
arranged in non-increasing order), and 
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Thus the eigenvalues of W° + B are shifted down relative to those of W° (the greatest eigenvalue 
of B is negative or zero), and similarly, the eigenvalues of W° + B + C — W K are shifted down 
relative to those of W° + B. Hence, 
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and, substituting into (19), 



£{b) < 0, 

the quantum mechanical effect must be to lower the kink energy. 
We have assumed, on physical grounds, that 



£ N {b) 



y2hrn ^~^ 



(fif(6)-n°), 



(28) 



(29) 



the Casimir energy of the order N truncated system converges to a constant value as N grows 
large. Numerical evidence for this is presented in figure 1, a graph of N against £jy(0)/m for various 
values of hm. Convergence is very fast for large hm (the hm =1.9 curve is essentially flat for 
N > 3) but much slower for small hm. Of course, this happens because the more finely space is 



discretized, the more degrees of freedom the kink structure is spread over. There is also a much 
smaller exacerbating effect due to the hm dependence in the kink solution, (ra) — large hm kinks 
are sharper in "real" x space than small hm kinks. Note that, as proved above, the Casimir energy 
is always negative and that its magnitude grows large for small hm as we expect from consideration 
of the unrenormalized continuum model M . Although 6 = was chosen for these data, the rates of 
convergence are virtually independent of b £ [— h/2, h/2]. 
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Figure 1: Convergence of the Casimir energy. 
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Figure 2: Position dependence of the Casimir energy. 



The Casimir energy of suitably truncated systems is plotted as a function of kink position b in 
figure 2. It attains its maxima where the kink is located exactly on a lattice site and its minima 
where the kink is halfway between two sites. For large hm the resemblance to a sinusoidal curve 
is remarkable: figure 3 shows a sinusoidal fit for hm = 1.8. The resemblence deteriorates as hm 



decreases because numerical errors become relatively large as the Casimir energy barrier [£ (0) — 
£(h/2)]/m decreases (figure 4). Striking though the similarity is, the author can find no analytic 
evidence to support the claim that £(b) really is sinusoidal. 
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Figure 3: Sinusoidal fit to the hm = 1.8 Casimir energy. 
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Figure 4-' The effect of discreteness on the Casimir energy barrier. 



It is the Casimir energy barrier (figure 4) which is physically most relevant and which will 
directly affect classical kink dynamics. In the continuum limit full continuous translation symmetry 
is recovered, so the barrier should disappear, as is indicated by the plot. In fact, in common with 
other (classical kinetic) discreteness effects of the LSGM [g, the barrier is very small for all hm < 1. 
If £(b)/m for a given b/h is a strictly increasing function of hm £ (0, 2), as it appears to be, then 
given that £ (b) < 0, the barrier should approach a finite value in the hm — ► 2 limit. Numerical 
evidence suggests that it does and that this value is around 0.17. 



5 Concluding Remarks 

We have seen that quantum fluctuations around the kink configuration spoil the level kink val- 
ley bottom by introducing a Casimir energy which depends periodically on the (classical) kink 
position. This energy has been computed numerically in the weak coupling approximation and 
found to be approximately sinusoidal, maximum when 6 = 0, ±/i, ±2/i, . . . and minimum when 
b = ±h/2, ±3/i/2, ±5/i/2, . . ., the difference between these extrema being large for large hm but 
rapidly vanishing in the continuum limit. It is superficially similar to the Peierls-Nabarro poten- 
tial of the Frenkel-Kontorova model || (the conventional lattice sine-Gordon model) but is entirely 
different in origin, being a purely quantum effect. Since the Casimir effect is a genuine physical 
phenomenon, experimentally verified in the context of Quantum Electrodynamics, we are led to the 
conclusion that classical kinks in this lattice model may be "pinned" by the quantum mechanics of 
the orthogonal modes. 

Finally, following a suggestion by Gibbons and Manton |J we could attempt to include the effect 
of variation of orthogonal mode frequencies on the quantized geodesic approximation of kink motion 
by including the Casimir energy as an extra potential term in the Hamiltonian: 

H GA = --A + S(b) (30) 

where A is the covariant Laplacian on the submanifold of static kinks. Given the periodic nature 
of £, we would then expect band structure in the kink spectrum. There are two objections to this. 
First, the kink kinetic terms are of order A 2 (due to the large kink mass), while the potential was 
expanded only up to order A , so the perturbative expansion is not consistent. Second, picking out 
the — ^A kinetic term at order A 2 is also not consistent because there are several other kinetic terms 



at this order neglected in equation (30). This may be seen by reexpressing the kinetic Hamiltonian 
of the full quantum field theory, 

h^ dtp 2 . 
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(an infinite dimensional Laplacian) in terms of the kink translation mode b and the orthogonal 
modes £j. The resulting formula is very messy, and includes order A 2 cross terms with derivatives 
d 2 /d^jdb and d 2 /d£jd£k in addition to the terms included in (|30|). The analogous expression in the 
continuum model || can be greatly simplified by boosting to the kink's rest frame, but this trick 
is not available to us here. These considerations cast doubt on the suggested procedure, at least in 
this case. However, one should note that a naive expansion of the Hamiltonian in A may not be 
the most physically relevant procedure when considering quantum kink dynamics. For example, one 
could imagine making the demand that the kink's "speed" (in an appropriate quantum sense) be of 
order A . Given the kink's large mass in the weak coupling approximation, one would then expect 
the kink kinetic term to dominate the Hamiltonian. So the legitimacy of the suggested procedure 
remains an open and somewhat controversial question. 
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